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The stability of color-flavor locked (CFL) strangelets is studied in the three-flavor Nambu-Jona- 
Lasinio model. We consider all quark flavors to be massless, for simplicity. By making use of the 
multiple reflection expansion, we explicitly take into account finite size effects and formulate the 
thermodynamic potential for CFL strangelets. We find that the CFL gap could be large enough 
so that the energy per baryon number of CFL strangelets is greatly affected. In addtion, if the 
quark-quark coupling constant is larger than a certain critical value, there is a possibility of finding 
absolutely stable CFL strangelets. 
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I. INTRODUCTION 



The properties of quark matter have attracted a good deal of interest since it has been suggested that strange 
quark matter could be a ground state of strongly interacting matter. 0.0 The most likely place for strange quark 
matter, perhaps, is the interior of compact stars even if it is not absolutely stable. The existence of absolutely stable 
strange quark matter is still an open question and it may be realized in the form of strangelets (small lumps of 
strange quark matter with roughly the same amount of up, down and strange quarks). Theoretically, there has been 
various investigations of the stability of (non)strange quark matter. Within the MIT bag model, Farhi and Jaffe Q 
found a reasonable window of the model parameters (i.e. the bag constant, current quark masses and strong coupling 
constant) for which strange quark matter is stable, while nonstrange quark matter is unstable as compared to a gas 
of 56 Fe. By including finite size effects (surface tension and/or curvature energy), similar analyses have been done for 
finite lumps of quark matter. 0,0, QUE El It is now well established that finite size effects increase the energy of finite 
quark lumps. 

During the last decade, significant advances have been made in our understandings of the phase structure of hot 
and/or dense QCD. 7] At the present time, it is widely accepted that a color-flavor locked (CFL) phase @ is the 
ground state of cold, dense quark matter. 1 Madsen 13] has studied the stability of CFL strangelets and found that 
CFL strangelets are significantly more stable than normal (unpaired) strangelets (see also Ref. 0]). However, that 
work ignored issues of density dependence of the bag constant and the CFL gap. What remains a question is the 
possible effects of the phase structure. In the previous paper |Tfij . we have studied the chiral and 2SC phase of finite 
quark lumps. We found that finite size effects enhance the restoration of chiral symmetry. In this case, it is likely 
that strangelets lie in the color superconducting phase. 

The purpose of this work is to study the behavior of the CFL gap in strangelets and to look at its effects on the 
stability of CFL strangelets. In order to describe the density dependent bag constant and the density dependent CFL 
gap we choose to use the three-flavor Nambu-Jona-Lasinio (NJL) model. (The NJL model is a simple tractable 
model to investigate the stability of quark matter/droplets [see, for example, Ref. H3]0 We consider up, down and 
strange quarks to be massless, for simplicity. In order to take account of finite size effects, we apply what is called 
multiple reflection expansion (MRE).[l8j The MRE has been used for calculating thermodynamic quantities such as 
the energy per baryon number and the free energy of finite quark lumps. As far as the general structure is concerned, 
the results are in good agreement with those given in the MIT bag model. 2, :G] Using the NJL model with the MRE, 
we formulate the thermodynamic potential for spherical CFL strangelets. We then investigate the stability of CFL 
strangelets, including the dynamical effects as well as the finite size effects. 

The organization of this paper is as follows. In the next section, we briefly review the thermodynamic potential of 
the CFL phase in bulk. Then, we formulate the thermodynamic potential for CFL strangelets. The gap equation and 
related thermodynamic quantities are also derived. In Sec. 3, we present numerical results. Finally, Sec. 4 is devoted 
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1 As discussed in Refs. I^ ll(llllill2i . it is of importance to take account of the constraint imposed by color and electric neutrality. For 
neutrino free, color and electrically neutral quark matter at zero temperature, it has been shown that the CFL phase opti mizes th e 
pairing energy and, then, is favored over the 2-flavor superconducting (2SC) phase in all (or almost all) the range of densities. Il& ITU 
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to the conclusions. 



II. THERMODYNAMIC POTENTIAL 



First, we briefly survey the thermodynamic potential of the CFL phase in bulk. To describe the CFL phase in the 
NJL model we follow Buballa and Oertel.[n| The U(3)l x U(3)/{ symmetric NJL Lagrangian is 

8 

C = qt-fdrt + Gi X! Wrf + &l5T l q) 2 ] , (1) 
z=0 

where q denotes a quark field with three flavors (Nf — 3) and three colors (N c — 3), the coupling constant G\ has a 
dimension [Gi] — [mass]~ 2 and the Gell-Mann matrices t % (i = 1, • • • ,8) with r° = ^/2/3J t act in the flavor space. 
In the case of the CFL, the quarks form the following (3 C , 3/) condensate 

Af = (q? gf) ex C-y 5 e^ x e ijX , (2) 

where i,j denote flavor indices, a, (3 denote color indices and repeated indices are summed. In principle, there must 
be the additional (6 C ,6/) condensate. However, we neglect such condensate because it is much smaller than (3 C ,3/) 
condensate ja l2(l | and its contribution to the thermodynamic potential is negligible. The interaction corresponding 
to the (3 C , 3/) CFL paring can be written as 

Ccfl=G 2 Yl (qil5T i X a Cf) {gCi^Wq) , (3) 

i=a=2,5,7 

where C is a charge conjugation matrix, defined by C _1 7 A1 C = — 7J and C T — — C, and t % (i = 2, 5, 7) [X a (a = 2, 5, 7)] 
denote the antisymmetric generators of SU(3)/ [SU(3) C ]. The coupling constant G2 can be obtained from Eq. 0} by 
making use of the Fierz transformation. However, we leave G2 a free parameter. In this approximation, all nine (three 
colors times three flavors) quarks participate in the pairing and yield eight quasiparticles with gap A and one with 
gap — 2A. (These quasiparticles correspond to an octet and a singlet of the unbroken SU(3), respectively.) Then, in 
the mean-field approximation, the thermodynamic potential fl = f2(A; /1, T) at finite quark chemical potential /x and 
temperature T is given by 
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- r E E f-^sH^+a), (4) 



± n— — 00 



where e± = yj{k± [i) 2 + A 2 , £± = y / (k± fi) 2 + 4A 2 and uj n = (2n + l)irT (n = 0, ±1, ±2, • • • ) denote the fermionic 
Matsubara frequencies. The summation over u n is straightforward. Henceforth, we restrict ourselves to T = 0. The 



T ->• limit of Q is 



± 



Here, we have introduced the ultraviolet cutoff Ayv m the three-dimensional momentum space. 

Now we incorporate finite size effects into the thermodynamic potential. To this end, we use the density of states 
derived from the MRE.0, 0, 0] In the MRE framework, the density of states for a spherical system is given by 
fc 2 PMRE/(27r 2 ), where 



Pmre = pMRE(k,m,R) 

kR Jb \mj (kR) 



6tt 2 , f k\ 12tt 2 , ( k ' 
1 + T77/S - + 777^2 fc -+•'•. ( 6 ) 



m 



with m being the Dirac mass of quark and R being the radius of the sphere. The second (third) term on the right-hand 
side represent the surface (curvature) contribution to the fermionic density of states. The ellipsis implies higher order 
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terms in 1/R, which are neglected throughout. In the case of massless fermions, the functions fs and fc have the 
following limits, 

lim fs(k/m) = , lim fc(k/m) = -l/(24^ 2 ). (7) 

m— *0 m— >0 

Then, we use the following MRE density of states, 

PMRE = 1 ~2(W' (8) 

One can see that the finite size effects reduce the density of states (the decreasing tendency is more pronounced at 
low momenta) and /?mre becomes negative at small kR. To avoid the unphysical negative density of states we shall 
introduce an infrared cutoff Air = \/2/ (2R) in the momentum space. 

Using the density of states ©, we express the effective potential of the spherical CFL strangelets as follows, 
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n M RE = -^-^/ (8e± + £t), (9) 



where we have introduced the following notation, 



MRE 



A UV k 2 dk 

TT^-PMRE- \W) 



in 



Here, we emphasize that the Fermi momentum fcp(= /i) in Eq. © is common to all nine quarks. In the case of the 
three flavors of massless quarks, a strangelet with the common Fermi momentum is automatically color and electrically 
neutral. 

A strangelet must be in a color singlet state. As noted in Refs. (lcLl2lj. color neutrality is nothing but a prerequisite 
condition for color singletness. However, it has been shown that, as long as a quark lump of color superconductor 
is color neutral, the projection over a color singlet state makes a negligible contribution to the thermodynamic 
potential. {21j Thus, color neutrality is a good approximation to color singletness. A small strangelet (A <C 10 7 ) 
does not need to satisfy electric neutrality. This is because the electron Compton wavelength is larger than such a 
strangelet and, then, electrons mainly stay outside of the quark phase. Q Therefore the phases of small strangelets 
could be different from that of quark matter in bulk. 

We focus on CFL strangelets embedded in a vacuum and derive a set of coupled equations with the MRE. For 
computation of a finite system, we choose a fixed radius R. First, we derive the gap equation, which is the extremum 
condition of ^mre with respect to A: 

«W 3A _4AW (± + Pl=0. (11) 



<9A 2G 2 ^7mre\ £ ± £± 

In addtion, we need to take account of the pressure balence relation. The relative pressure inside strangelet Pmre 
is given by 

-Pmre = Pmre — Pvac, (12) 

where Pmre = -Omre and p vac has been introduced to measure the pressure relative to outside of strangelet. Note 
that Pvac is nothing but the pressure of chirally broken vacuum. In the mean-field approximation, its value can be 
determined by using the Lagrangian as follows: 

/• A uv (pk i 3M 2 

Here, M is the dynamically generated quark mass, which is a nontrivial solution to the following equation, 

3M —r\ d \ = =o d4) 

4(17! f C J (2tt) 3 Vfc 2 + M 2 V ' 

We solve Eqs. (|ll|l and i|12fl self-consistently and, then, compute the baryon number of the strangelet A by making 
use of the following relation, 



A = Vn B , 



(15) 



4 



> 
0) 



0.11 



0.1 



0.09 



0.08 



0.07 




0.3 



0.6 



FIG. 1: The CFL gap A for finite volumes. The solid curves are obtained for three values of R: R = 3, 5 fm (from bottom to 
top). The dotted line refers to the case of R — > oo. We have taken G2/G1 = 3/4. 



where V — 4ttR 3 /3 is the volume of the spherical strangelet and 
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(16) 



is the baryon number density of the strangelet. 



III. NUMERICAL RESULTS 

In numerical calculation, we use the following set of parameters: Atjv = 0.6 GeV, G1A 2 = 2.311. [19j At T = fi = 0, 
these parameters yield the constituent mass M = 0.35 GeV and the bag constant B 1 / 4 = 0.182 GeV. The quark-quark 
coupling constant G2 has been set to be G2/G1 = 3/4 in the literature. 0, 0, |2^ This relation between G\ and G2 
is obtained from a four-fermion interaction with a quantum number of a one-gluon exchange, 



£oGE = -9^(g7 A1 A Q g) 2 , 



(17) 



by making use of the Firez transformation. Moreover, this is consistent with that determined by fitting the nucleon 
mass within a Fadeev approach. j2^ However, we vary G2 in the range of G2/G1 = (1/2 — 1) to see its effect on the 
stability of strangelets. 

We begin with the discussion of the solution to the gap equation (|llf> . Figure 1 shows the CFL gap for finite 
volumes as a function of the quark chemical potential. We find that the finite size effects reduce the gap. However, 
the decrease in the gap is not pronounced even at small radii; therefore, the R — > 00 limit is a good approximation as 
to the size of the gap. 2 This is easily understood as follows. Color superconductivity is brought about by the quarks 



2 In this respect the result is similar to that in Ref. |2 ll where finite size effects on the 2SC gap have been studied by confining the quarks 
to a cubic box. 
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FIG. 2: The CFL gap A as functions of the baryon number A. The curves correspond to the cases of G2/G1 = 1/2, 3/4, 1 
(from bottom to top). 



near the Fermi surface. Recall that at a fixed radius R the decrease in the density of states is not pronounced at large 
momenta. Then, if the Fermi momentum is sufficiently large, the finite size effects do not greatly affect the gap. 

In order to look at CFL strangelets, we solve the gap equation the pressure balance relation l|12ll and the 

baryon number relation 1151) . self-consistently. By numerically solving these equations, we obtain the baryon number 
dependence of the gap, the quark chemical potential and the radius. After that the energy per baryon number of the 
pressure-balanced strangelet, 



E 
A 



p=0 



£_ 

n B 



3/z, 



(18) 



p=0 



can be evaluated. In Eq. (|18|) . £ — JImre + ^J2 a na ~ £vac denotes the energy density of the strangelet, where 
a G (u,d,s), n u — rid — n s (~ tib) are the quark number densities, and e vac (= — p V ac) is the energy density of the 
chirally broken vacuum. 

Figure 2 shows the CFL gap as a function of the baryon number for the cases of G2/G1 = 1/2, 3/4, 1. Each curve 
remains approximately constant as long as A is not too small. The change of the gap at small baryon numbers 
(A ;$ 100) is a consequence of the fact that the actual quark chemical potential increases when A decreases. We note 
that the results for A ;$ 10 should not be considered robust ones, because in this regime the quark chemical potentials 
are close to the ultraviolet cutoff Atjv- 

In Fig. 3, the energy per baryon number E/A for CFL strangelets is plotted as a function of the baryon number. 
For comparison, we also present the result for normal strangelets. As is evident from the figure, CFL strangelets 
are considerably more stable than normal strangelets. Needless to say, this decrease in E/A is due to the pairing 
energy contribution to the thermodynamic potential Omre- Note that E/A depends on the size of the gap; hence, 
the quark-quark coupling constant G2. As G2 is increased, the gap increases and then the contribution to JImre from 
the pairing energy also increases. Thus, it is reasonable that the growth of G2 tends to lower E/A of CFL strangelets. 
The gap can be large enough to have a great effect on E/A. For instance, we observe that at G2 = G\ CFL strangelets 
of A > 160 are absolutely stable. Within our model parameters (Atjv and Gi), G2/G1 needs to exceed 0.88 for the 
existence of the absolutely stable CFL strangelets. Note also that each curve shows the typical behavior of E/A of 
finite quark lumps. E/A(= 3/Lt) increases with decreasing A. (As mentioned earlier, this is a consequence of the 
increase of the quark chemical potential by the finite size effects.) On the other hand, as A grows, E/A approaches 
the value without the finite size effects. 
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FIG. 3: The energy per baryon E/A for CFL strangelets. The solid curves correspond to the cases of G2/G1 = 1/2,3/4,1 
(from top to bottom). The dotted line refers to E/A for normal strangelets. 



In this work, we have explored the color-flavor locked phase in strangelets. We used the three-flavor Nambu-Jona- 
Lasinio model. To take account of finite size effects we applied the multiple reflection expansion. Using the NJL model 
with the MRE, we formulated the thermodynamic potential for CFL strangelets and derived related thermodynamic 
quantities. We clarified the behavior of the CFL gap and its effects on the stability of strangelets, although our model 
is nothing but a toy model, neglecting the finite strange quark mass. The CFL gap is almost independent of A as 
long as A is not too small. Due to the contribution from the formation of the CFL condensate, CFL strangelets are 
more stable than strangelets without the CFL. Further, if G2 is larger than the critical value, we have a chance to 
find absolutely stable CFL strangerets. These results complement some of the conclusions made by Madsen.^l Of 
course we cannot completely accept the existence of absolutely stable CFL strangelets, because we do not have a 
good knowledge of Gi. However, a Fadeev approach to baryons [23l |24| would be a guide to the determination of 
G2. Moreover, we had neglected many aspects of dense QCD. Most important is probably the finite strange quark 
mass m s . When m s is nonzero, we have to take account of unlocking phase transition. [25j As pointed out in Ref. 
fl2| . the main reason for the difficulty to find the stable 2SC matter arises from the constraint of electric neutrality. 
In contrast, small strangelets, which can be regarded as free from this constraint, could be in the 2SC phase. |26| It 
should be mentioned, however, that the increase of the chemical potential caused by the finite size effects might prefer 
the CFL phase to the 2SC phase. It would be of great importance to study competition with oth er p hases (hadronic 
phase, 2SC phase, and so on), taking account of nonzero m s as well as color neutrality. |9UlOUllUl2| A more careful 
analysis is left to future studies. 

It should be also noted that the MRE contains several problems concerning its reliability. First, the density of 
state pmre should have the terms proportional to l/R 3 , 1/-R 4 , and so on. These terms would be dominant at small 
radii, i.e., at small baryon numbers. Furthermore, pmke in Eq. 10 causes the unphysical negative density of states 
at small kR. Although our qualitative results for relatively large baryon numbers (A 2; 100) would not change, it is 
desirable to take more rigorous way of including finite size effects. In particular, as noted in Ref. 0, a detailed study 
of strangelet decay modes will have to rely on shell model calculations. 

In this work, we restricted ourselves to zero temperature. However, for applications such as heavy ion collision 
experiments and the cosmological quark-hadron phase transition, it is interesting to investigate strangelets at finite 
temperature. The present approach is straightforwardly generalized to nonzero temperature. Studies along this line 



IV. CONCLUSIONS 
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are now in progress. 
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